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payload  was  launched  from  Wallops  Is.,  Virginia  and  the  release 
occurred  at  348  km  altitude  at  5:20  UT  on  9  Nov  1983. 

Subsequent  sections  of  the  report  discuss  methods  used  for 
determination  of  video  camera  pointing  from  the  recorded  star 
fields,  reduction  of  range  to  the  glow,  analysis  to  determine 
growth  of  the  cloud,  and  cloud  motion  (atmospheric  wind) 
derivation  from  its  time  dependent  positions. 

2.  Determination  of  Camera  Orientation 

Triangulation  camera  orientation  is  normally  determined 
very  accurately  using  vector  mathematics  to  relate  image  plane 
locations  of  star  images  with  their  known  equatorial  coordinate 
system  positions  (right  ascension  a,  and  declination  6). 

Data  for  several  stars  and  their  images  are  processed  utilizing 
the  methods  of  Ref's  2-3  to  establish  the  focal  length  of  the 
camera  and  the  direction  of  the  optic  axis  in  space.  The 
direction  vector  is  intlally  calculated  in  the  equatorial 
system  but  is  easily  expressed  in  other  coordinate  systems  by 
means  of  coordinate  transformations  (rotations  about  coordinate 
axes) .  Typical  precisions  (standard  deviations)  of  the  measure¬ 
ments  for  a  short  focal  length  lens  (f  -  25  mm)  are  0.05  cm  for 
focal  length,  and  0.01°  for  azimuth,  elevation,  and  horizonal 
tilt  (orientation  expressed  in  the  horizon  coordinate  system 
where  the  horizonal  tilt  is  the  angle  between  the  camera  x-axis 
and  the  horizon).  Note  that  procedures  which  are  normally 
applied  when  lens  distortion  is  negligible  (i.e.  a  single 
focal  length  completely  describes  the  image  field)  were  ex¬ 
panded  and  adapted  for  the  90°  lens  used.  Furthermore,  the 
small  number  of  stars  visible  reduced  the  precision  of  the 
camera  orientation. 

Digitization  of  the  video  and  determination  of  image  plane 
coordinates  of  star  and  glow  images  was  performed  on  the  AFGL 


Figure  1.  Video  image  of  the  7774A  glow  at  0520  UT, 
9  Nov  1983  (40  sec  after  release).  Note 
that  many  of  the  dot  Images  are  not  stars 
(also  in  Fig's  2-3). 


•/ 

Figure  3.  Video  image  of  the  7774A  glow  at  0522  UT, 
9  Nov  1983  (160  sec  after  release). 
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While  such  large  errors  are  usually  indicative  of  one  or 
more  stars  which  have  been  misidentified,  the  problem  here 
appears  to  be  associated  with  a  (presumably  radially  symmetric) 
off-axis  distortion  commonly  found  in  very  wide  angle  camera 
systems  associated  with  its  variation  in  focal  length  with 
angular  distance  from  the  optic  axis.  Compensation  for  this 
lens  distortion  was  attempted  several  times  with  little  or  no 
improvement  in  the  resulting  optic  axis  direction  vector  error. 
The  distortion  curve  used,  shown  in  Fig  4,  was  taken  from  a 
90°  field,  0.345  cm  focal  length  lens  of  broadly  similar  design 
and  shifted  to  coincide  with  the  measurable  distortion  at  the 
locations  of  identifiable  stars.  Stars  used  to  measure  this 
distortion  were  digitized  from  a  video  frame  taken  through  a 
6300A°  filter  where  many  more  stars  are  identifiable  than 
in  the  the  near-lR  7774°  images.  (There  were  still  insufficient 
stars  to  adequately  determine  the  distortion  curve.)  Note 
that  the  optical  system  used  with  the  image  lntenslfler  camera 
has  a  relatively  complex  telecentric  lens  similar  to  Fig  5,  with 
the  complexity  deriviving  mainly  from  the  small  ray  divergence 
needed  at  the  position  of  the  interference  filter  to  obtain  a 
nearly  monochromatic  response  at  all  off-axis  angles. 

A  further  factor  that  affected  the  distortion  correction 
was  Inability  to  accurately  locate  the  center  of  the  circular 
image  (optic  axis).  That  is,  the  correction  as  applied  may 
have  been  asymmetric.  (The  optical  system  lacked  fiducial 
markers  for  establishing  a  coordinate  system  in  the  film  plane, 
and  hence  the  optic  axis) .  We  located  the  center  by  reference 
to  several  points  on  the  (relatively  diffuse)  edge  of  the 
circular  field  of  view.  Primarily  the  vertical  and  horizontal 
diameters  were  bisected  to  obtain  the  center  point,  since 
obvious  irregularities  exist  over  most  of  the  circumference. 


Mislocation  of  Che  optic  axis  not  only  affects  the  distortion 
correction,  but  also  shifts  coordinate  locations  (relative 
to  the  axis)  and  adversely  affects  coordinate  transformations 
(rotations  about  coordinate  axes)  performed  during  orientation 
and  triangulation  calculations. 

In  order  to  establish  the  optic  axis  direction  vector  in 
a  manner  which,  while  not  precise,  would  be  consistent  from 
frame  to  frame,  it  was  decided  to  use  only  the  two  stars 
nearest  to  the  center  in  each  case.  The  best  estimate  of  the 
image  plane  position  of  the  optic  axis  was  used  and  the  assump¬ 
tion  made  that  there  was  zero  lens  distortion.  The  star  image 
digitization  program  which  predicts  the  image  plane  location 
of  stars  in  the  field  of  view  was  then  run  for  several  focal 
lengths  (0.28  <  f  <0.34  cm).  The  focal  length  resulting 
in  the  best  match  (f  -  0.31  cm)  of  predicted  with  actual 
measured  positions  for  the  two  star  images  was  chosen  as  the 
effective  focal  length  for  use  in  determining  the  direction 
vector. 

Finally,  the  direction  vector  was  determined  for  each  frame 
using  these  same  two  stars.  These  are  shown  in  Table  I 
expressed  in  the  horizon  system.  Although  no  standard  devia¬ 
tions  may  be  calculated  from  a  two  star  determination  some 
worst  case  (assuming  ±  1  pixel  error  in  star  coordinates) 
estimates  of  error  for  these  data  are  also  shown.  More  real¬ 
istic  estimates  based  upon  the  interpolation  used  to  locate 
the  star  image  coordinates  are  on  the  last  line  of  the  table. 
(Interpolations  were  done  to  1/10  pixel  and  a  ±  1/5  pixel 
error  assumed) .  The  angular  errors  translate  to  spatial 
errors  of  '1.2  km  at  ranges  of  -350  km.  Since  the 
release  image  is  located  near  the  center  of  the  frame,  the 
large  tilt  (rotational)  error  produces  only  a  very  minor 
error  in  spatial  location  of  the  glow. 
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t.  Triangulation  and  Velocity 

Since  images  of  the  release  were  only  available  from  one 
observation  site,  a  single  site  triangulation  (projection) 
was  used  to  determine  the  spatial  position  of  the  glow  as  a 
function  of  time.  This  method  assumes  that  the  altitude  of 
the  glow  is  known,  scales  the  line  of  sight  vector  to  the  glow 
to  intersect  the  assumed  altitude,  and  calculates  the  latitude 
and  longitude  of  this  point. 

Briefly,  the  image  plane  glow  coordinates  were  converted 
to  a  vector  in  the  camera  system.  Then  utilizing  the  camera 
azimuth,  elevation,  and  tilt  (calculated  above  from  the  star 
field),  the  camera  location  (latitude,  longitude,  altitude), 
and  the  exact  time  that  the  image  was  exposed,  a  series  of 
coordinate  transformations  converted  the  camera  vector  to  a 


geocentric  coordinate  system  vector  L.  The  vector  pc  from 
earth  center  to  the  camera  was  then  calculated  from  an  ellip¬ 
soidal  earth  model.  An  iterative  procedure  then  added  Pc  + 

AL  -  pc  where  A  is  a  scale  factor  to  change  the  length  of  L 
and  pQ  is  the  vector  from  earth  center  to  the  glow.  As  each 
PG  was  obtained,  the  earth  model  was  used  to  calculate  the 
latitude,  longitude,  and  altitude  Z  defined  by  pq.  If  Z  did 
not  equal  the  assumed  glow  altitude  H,  A  was  modified  and 
another  Iteration  performed  until  Z  ■  H.  Latitudes  and 
longitudes  thus  determined  for  a  series  of  time  separated 
images  provide  the  basis  for  determining  glow  velocity. 
Altitude  precision  is  estimated  as  approximately  t  2  km 
based  upon  camera  orientation  errors. 

Table  II  shows  spatial  positions  determined  from  the 
images  based  upon  two  assumptions  concerning  the  vertical 
velocity  Vv  of  the  glow.  Included  in  the  data  is  the  rocket 
trajectory  position  at  the  time  of  release.  The  first  set  of 
data  assumes  that  Vv  “  0,  that  is  the  altitude  does  not  change 


with  time.  It  has  been  observed  however ,  that  the  altitude 
of  the  region  of  maximum  radiating  material  density  decreases 
with  time.  Computer  assisted  calculations  (Ref  1)  show  the 
fall  rate  to  average  '•112  m/sec  in  the  first  40  sec,  -46  m/sec 
in  the  next  60  sec,  and  -38  m/sec  in  the  following  60  sec. 

The  altitudes  used  for  the  triangulation  position  determination 
for  the  second  set  of  data  in  Table  II  were  obtained  from  these 
calculations  and  are  also  shown  in  Fig  6.  Horizontal  velocities 
were  calculated  from  the  position  differences  and  are  shown  for 
each  time  interval.  Average  velocity  over  the  total  interval 
has  been  determined  from  least  squares  analysis  of  positions. 
Position  north  and  south  relative  to  release  point  and  the 
velocity  for  each  time  interval  are  shown  in  Fig's  7  and  8. 


Figure  7.  Horizontal  position  of  the  glow  relative  to  release  point 
at  times  0520,  0521,  0522  UT  assuming  that  the  glow  region 
is  decreasing  in  altitude. 


4.  Growth  of  the  SFfi  Induced  Glow 

The  rate  of  growth  of  the  image  of  a  tracer  glow  may  be 
used  to  determine  the  rate  of  diffusion  of  the  tracer  (or  its 
reaction  product (s))  into  the  ambient  atmosphere  (Ref  6).  Only 
relative  exposure  (irradiance  x  time)  need  be  recorded  in  the 
image  (photographic  or  electronic)  to  determine  the  diffusion 
coefficient.  Changes  in  the  "gain”  of  the  recording  system 
(lens  aperture  or  electronic  amplification)  may  also  effect  the 
gaussian  width,  in  that  the  slope  of  the  curve  of  display  in¬ 
tensity  vs  faceplate  irradiance  would  change.  The  assumption 
was  made  that  the  response  of  the  image  intensifies  and 
vidicon  to  irradiances  was  linear  over  the  small  brightness 
range  exhibited  at  7774a°. 

The  recorded  image  is  a  composite  of  exposure  due  to  the 
glow  and  that  due  to  background  (night  sky  in  this  case). 

This  background  must  be  subtracted  leaving  only  the  glow 
exposure  itself,  which  is  typically  gaussian,  to  be  analyzed. 
The  diffusion  coefficient  is 

D  -  (h2  -  hQ2)/4(t  -  tQ) 


where  h  and  hQ  are  the  glow  exposure  halfwidths  at  time  t  and 
tQ.  For  a  gaussian  glow  h2  is  best  determined  from  a  least 
squares  fit  of  An  E  and  r2  where  E  is  the  glow  exposure 
at  distance  r  from  the  glow  center.  The  slope  derived  from 
this  fit  is  -1/h2. 

The  equation  for  D  has  been  used  extensively  even  when  the 
glow  is  not  strictly  gaussian  (the  least  squares  analysis 
straight  line  will  not  necessarily  be  indicative  of  h2), 
however  h  is  measured  directly  from  the  glow  image  distribution 
at  the  point  where  exposure  ■  1/eth  of  the  peak  glow  exposure. 
For  a  gaussian  distribution  there  is  no  difference  between 
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Table  III.  Measured  and  averaged  values  of  SFg  glow  halfwidths. 

An  average  slant  range  of  336  km  and  focal  length 
of  0.31  cm  were  used  to  obtain  spatial  halfwidths. 


TIME 

(sec) 

42.5 


102.5 


162.5 


42.5 


102.5 


162.5 


h  IN  IMAGE  h  AT  GLOW 
PLANE  (cm)  (cm) 


.01835 


.02802 


.03695 


.01618 


.02474 


.02753 


1.98x10* 

3.04x10* 


4.00x10* 


1.75x10* 

2.68x10* 


2.98x10* 


h2  AT  GLOW 
(cm2) 

3.95x1012 

HORIZONTALLY  IN 
9.28x1012  IMAGE  PLANE 

1.60x1013 


3.07x1012 

VERTICALLY  IN 
7.19x1012  IMAGE  PLANE 


8.90x1012 


42.5 


1.86x10* 


3.47x1012 


102.5 


2.86x10* 


8.17x1012  AVERAGE 


162.5 


3.49x10* 


1.21x1012 
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5.  Conclusions  and  Recommendations 

Image  Intensified  video  records  of  the  7774a  01  emission 
Induced  by  a  release  of  sulfur  hexafluoride  molecules  Into 
the  atmosphere  have  been  analyzed  using  the  AFGL  Coordinate 
Digit Izer/lmage  Analyzer.  Camera  axis  direction  vectors  were 
obtained  by  analyzing  the  recorded  star  fields  and  spatial 
position  of  the  glow  and  its  range  from  the  camera  determined 
by  single  site  triangulation  (projection).  The  atmospheric 
wind  (velocity  of  the  glow)  was  deduced  from  a  time  sequence 
of  cloud  positions.  Further  analysis  was  conducted  to  determine 
the  growth  rate  of  the  released  material  into  the  ambient 
atmosphere,  from  which  an  effective  diffusion  coefficient  was 
derived.  No  geomagnetic  confinement  of  the  glow  species  was 
detectable.  With  absolute  radiometric  calibration  data,  we 
would  be  able  to  assess  the  absolute  surface  brightnesses  of 
the  glow  (from  which  column  densities  of  reactants  could  be 
inferred.) 

Use  of  this  intensified  and  filtered  imaging  system  has 
provided  the  first  measurements  of  atmospheric  wind  and 
diffusion  at  thermospheric  altitudes  above  '250  km.  In 
addition  it  provided  information  on  the  atmospheric  chemical 
reactions  of  SFg  and  its  successor  species.  Two  such  cameras, 
widely  separated,  would  allow  true  triangulation  to  such  glows 
and  more  precise  wind  measurements.  A  measurement  of  the  lens 
focal  length  as  a  function  of  off-axis  angle  would  improve  the 
triangulation  accuracy,  and  application  of  the  absolute  photo¬ 
metric  response  of  the  camera  would  result  in  still  further 
useful  information  about  both  the  ambient  atmosphere  and  the 
reactions  resulting  from  the  chemical  release. 
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